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ABSTRACT 

Quantitative dynamical models of galaxies require deprojecting the observed sur- 
face brightness to determine the luminosity density of the galaxy. Existing deprojection 
methods for axisymmetric galaxies assume that a unique deprojection exists for any 
given inclination, even though the projected density is known to be degenerate to the 
addition of "konus densities" that are invisible in projection. We develop a deprojection 
method based on linear regularization that can explore the range of luminosity densities 
statistically consistent with an observed surface brightness distribution. The luminosity 
density is poorly constrained at modest inclinations (i ^ 30° ) , even in the limit of van- 
ishing observational errors. In constant mass-to-light ratio, axisymmetric, two-integral 
dynamical models, the uncertainties in the luminosity density result in large uncertain- 
ties in the meridional plane velocities. However, the projected line-of-sight velocities 
show variations comparable to current typical observational uncertainties. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 
1439, NGC 7619) — galaxies: kinematics and dynamics — galaxies: structure 

1. INTRODUCTION 

Models of elliptical galaxies seek to understand the spatial distribution of their stars (e.g. 
Franx, Illingworth, & de Zeeuw 1991; Statler 1995), the structure of their orbits (Dehnen & Gerhard 
1993, 1994; Arnold, Robijn, & de Zeeuw 1995), and to secure evidence for the presence of dark 
halos (Saglia et al. 1993; Carollo et al. 1995) or black holes (van der Marel et al. 1994; Dehnen 

1995) . These problems can be probed using the observed velocities of the stars (van der Marel 
1991; Carollo & Danziger 1994), X-ray emission (Fabbiano 1989, 1995), gaseous disks and rings 
(Schweizer, Whitmore, & Rubin 1983; Whitmore, McElroy, & Schweizer 1987; Bertola et al. 1991; 
Franx, van Gorkom, &: de Zeeuw 1994), and gravitational lenses (Maoz &: Rix 1993; Kochanek 1995, 

1996) . Most modern dynamical models employ axisymmetric, two-integral models in which the 
distribution function depends only on the energy and the angular momentum about the symmetry 
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axis (e.g. Binney, Davies, & Illingworth 1990; Hunter & Qian 1993). Such models have simple 
solutions to the Jeans equations (Satoh 1980; see also Binney & Tremaine 1987 §4.2), and in some 
cases the distribution function can be determined (Dehnen & Gerhard 1993, 1994; Qian et al. 1995). 
Quantitative application of the two-integral models to a real galaxy relies on the deprojection of 
the observed surface brightness of the galaxy to determine its three-dimensional luminosity density. 

The uniqueness of the deprojection of a galaxy depends on the symmetries of the density. 
Spherical galaxies always have unique deprojections, and ellipsoidal galaxies with fixed axes have 
unique deprojections for known inclinations (e.g. Binney 1985). General axisymmetric galaxies 
do not have unique deprojections, except when the symmetry axis is in the plane of the sky 
(inclination i = 90°, Rybicki 1987). The projection operation destroys information about the 
Fourier components of the density which lie in a "cone of ignorance" of opening angle 90° — i about 
the symmetry axis. A "konus density" (Gerhard & Binney 1996), whose Fourier transform is non- 
zero only inside the cone of ignorance, can be added to the luminosity density without changing the 
projected surface brightness. Kochanek & Rybicki (1996) developed methods to produce families 
of konus densities with arbitrary equatorial density distributions. Simple konus densities generally 
look like "disks" because of the conical symmetry of the model in Fourier space; this is consistent 
with Rix & White's (1990) observation that disks would be nearly invisible in ellipticals, unless 
close to edge-on. 

We need to explore two issues to understand the effects of the deprojection degeneracy on 
axisymmetric models of ellipticals. The first issue is the existence of smooth konus densities. 
A real galaxy must have a smooth, monotonic, positive definite density profile. These physical 
restrictions limit the allowed range of konus densities. For example, none of the analytic konus 
densities of Gerhard & Binney (1996) or Kochanek & Rybicki (1996) are physical because they 
have unacceptable angular profiles at large radii due to the conditions imposed to find analytic 
solutions. The second issue is that even if the luminosity density is underconstrained, the degeneracy 
is dynamically interesting only if the inferred velocity dispersions or velocity dispersion profiles are 
uncertain by more than the observational errors. 

The projection operator P is a linear operator I = Pzv between the luminosity density v and 
the surface brightness I. Deprojection corresponds to inverting the linear operator v = P -1 /. 
Methods to invert P must cope with three generic problems: true degeneracies, limited sampling, 
and amplification of noise. True degeneracies, such as the konus densities, correspond to eigenvalues 
of P that are zero, and they can be handled only by adding additional constraints on the inver- 
sion such as smoothness, positivity, or monotonicity. The surface brightness is always discretely 
sampled, but we want to determine the continuous luminosity density. If we try to determine the 
luminosity density at more points than are sampled in the surface brightness, the luminosity density 
is underconstrained by the data and additional constraints are required to perform the inversion. 
Most inversions are also ill-conditioned, because some of the eigenvalues of P are nearly zero. Small 
eigenvalues of P are large eigenvalues of P _1 , and they amplify small fluctuations of / into large 
fluctuations in v. The inversion method must suppress these oscillations. 
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Previous approaches to deprojection produce "unique" inversions because they have additional, 
hidden constraints. Lucy's (1974) method is the standard deprojection algorithm used in stellar 
dynamical models (e.g. Newton & Binney 1987; Binney et al. 1990; van der Marel, Binney, & 
Davies 1990; Gerhard 1991; van der Marel 1991; Dehnen 1995). It is a simple, iterative scheme 
that converges toward a density distribution that exactly fits the data; a "unique" solution results 
because of an implicit, nonparametric bias. Lucy's method also introduces numerical instabilities 
into the solution, so the iterations are manually halted at some point when the density still "looks 
good". Bendinelli's (1991) method expands the surface brightness in terms of Gaussian profiles, 
and then numerically fits the density to these terms. Palmer's (1994) method finds a finite angular 
polynomial series for the density whose projection fits the observed surface brightness to arbitrary 
accuracy. Functional fitting methods are fast, and frequently easy to program, but they are limited 
to a subset of the possible solutions. 

None of these existing methods provides a way of studying the true degeneracies of deprojection 
and their effects on dynamical models of axisymmetric galaxies. In §2, we develop a deprojection 
method based on linear regularization that allows us to explore the degeneracies of the inversion 
while keeping the density well-defined and physical. In §3, we test the algorithm on artificial galaxy 
images. In §4, we give results for two galaxies (NGC 1439 and NGC 7619), and in §5 we summarize 
our conclusions. Two appendices add some details of the numerical algorithms. 



2. METHODS 

2.1. Projection Geometry and Algorithm 

An axisymmetric galaxy has density distribution v(X, Y, Z), where the Z-axis is the symmetry 
axis of the galaxy. The surface brightness is I(x,y), where (x, y) are the coordinates in the plane of 
the sky, and the z-axis lies along the line of sight. The X- and x-axes coincide, and correspond to 
the line of nodes at the intersection of the plane of the sky and the equatorial plane of the galaxy 
(thus, for an oblate axisymmetric galaxy, the x-axis corresponds to the observed major axis of 
the galaxy). Because the galaxy is axisymmetric, its density distribution is completely specified by 
v(R, Z), where R = \J X 2 + Y 2 . For simplicity, we assume that the galaxy has a reflection symmetry 
about the equatorial plane, v{Z) = v{—Z). The surface density is related to the luminosity density 
by the projection operation, 

/+oo 
u(R,Z)dz. (1) 
-oo 

The transformation between the galaxy's coordinates and the observer's coordinates is 

X = x 

Y = y cos i — z sin i 
Z = y sin i + z cos i 
r 2 = R 2 + Z 2 = zu 2 + z 2 , (2) 
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where m 2 = x 2 + y 2 , and the inclination i = 90° corresponds to an "edge-on" galaxy with no cone 
of ignorance. 

The density of the galaxy v{R, Z) is divided into zones, Uj k = u(9j,rk), where one quadrant of 
the meridional plane {R, Z) is divided into N r radial zones and into iV a angular zones. The radial 
zones are logarithmically spaced in r from inside tz7 min to outside TX7 max , where m m i n and m max 
are the radii of the inner and outer surface density measurements. The angle 9 is the standard 
spherical polar angle with 6 = on the symmetry axis; the angular zones are equally spaced in 9 
(9j = j ■ AO, where j = [0 • • • N a — 1] and A9 = ir/2N a ). Reflection symmetry about the equatorial 
plane is implicit in the model, although the symmetry could be trivially removed. 

Although the zones Vjk represent the real density, the projection algorithm considers the galaxy 
to be made of "stacked blocks" in order to force the density to vary monotonically in anglefil. For 
an oblate galaxy, the total density of each successive zone increases with its angle 0j, and each 
additional contribution can be thought of as a block of density Vj k = p 2 k stacked on top of the 
previous blocks, with total density 



i=j i=j 

vjk = J2vik = J2p 2 i k - (3) 

i=0 i=0 



By construction, a density model expressed in terms of pjk is positive definite, oblatea, and mono- 
tonic in angle. 

The projected intensity is sampled at discrete points I\ m = I(xi m , yi m ), where I, m are arbitrary 
indices (e.g. corresponding to polar coordinates mi, ip m ). The contribution to ij m from density 
block i>jk is found by integrating along the line of sight through the block 

rz-2 

dljklm = / v jk (z)dz, (4) 

J Z\ 

where (z\,rx), (z2> r 2) ar e the coordinates of the intersection points of the line of sight from 
(xim,yim) with the edges of zone {jk}, and 

lira = X! dljklm- (5) 

jk 

We use a first-order projection scheme that linearly interpolates the density between radially- 
adjacent zones with 

j>jk(z) = Vjk + — —[rk-r(z)], (6) 
r k - r k+ i 



X A general two-dimensional density distribution that is monotonic in both angle and radius cannot be constructed 
in this way. 



2 For a prolate galaxy, the density decreases with angle Qy. Vjk = Xyj— 3 1 ^ik- 
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giving the first-order approximation 

Z2 ~ Zl 



dl. 



jklm 



z 2 r 2 - Z\r\ + w In 



Z2 + r 2 

zi + n 



jklm 



(7) 



(see Appendix A). Beyond the outermost density zones, the density is assumed to decrease as a 
power law, normalized to the density of the outermost zone at a given angle 9y. v oc z^o( r2 + s b)~ ab ^ 2 - 
The (small) contribution to the projection from this density "tail" is found by numerical integration 
of equation (4). 



2.2. Smoothing and Regular izat ion 

Given the projection, ij m , of the current model density distribution, and the observed surface 
brightness, we define a x 2 statistic for the goodness of fit, 

2 _ ( ~ 



Im 



&lr, 



x 2 = EP 2 — ^ . ( 8 ) 



where ai m is the noise associated with the measurement I* (we neglect any correlation function 
of the noise). If there are iV"d ata points I\ m , then a good fit should have x 2 — A'data with a one 
standard deviation error of Ax 2 = X 2 ~ ^data — i\/2^Vdata (f° r -^data ^> 1) if we neglect the number 
of degrees of freedom in the source model. The algorithm used to minimize the function requires 
the first and second derivatives of the x 2 with respect to the density kernels pi m ; these derivatives 
are given in Appendix B. 

Linear regularization is used to combat the problems that occur in unregularized solutions of 
integral equations such as equation (1). The primary problem is the tendency for the density to 
oscillate between the radial sampling points when there are more model density points than surface 
sampling points. The simplest solution is to find a good fit to the data while simultaneously 
minimizing a smoothing function 

hjk+l^jk+l — hjk v jk\ ^ 




hjk^jk J 

defined by the fractional variation in the radial density profile, weighted by the bias function 
hjk = h(rie, 9j). The bias function is used to weight the smoothness equally over all points, and to 
provide a bias slope at large radii where sky-level uncertainties poorly constrain the radial profile. 
The density should also vary smoothly with angle, so we add a second smoothing function 

H2 = ^{ ) ' <10> 

to prevent unrealistically sharp angular density variations. Ideally, h(r, 9) would be given self- 
consistently by the deprojected radial profile z/ -1 (r, 9). In practice, we fit a simple analytic model, 
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e.g. the power-law model 

I{w)=h{sl + w 2 f-^l\ (11) 

to the surface data points along the major axis, and analytically deproject it to give h(r); then we 
assume an constant axis ratio = b/a to give the angular variation h(rk,0j), e.g., 

hjk = 

Such radial surface brightness models qualitatively cover much of the observed range of galactic 
morphologies (Binney & Tremaine 1987). In addition to enforcing smoothness, the function h(r, 6) 
allows us to explicitly bias the profile to some arbitrary morphology. 

We then minimize the function F = x 2 + \H, where H = H± + nH2- The weighting factor 
k is somewhat arbitrary, but must be small enough to keep the radial profile acceptably smooth. 
For each image /* and bias function h, k is set by trial-and-error to be as large as is possible 
before significant radial fluctuations appear in the solution. The Lagrange multiplier A must be 
adjusted so that when F is minimized, the x 2 1S found to have a value in the range A^ata i V 2 A^ata ; 
where Adata is the number of surface density sampling points (we used solutions with \ 2 in t ne 
range Adata ± v^data)- The Lagrangian multiplier A is found iteratively: an initial A is chosen, the 
function is minimized, and the resulting x 2 value is linearly interpolated to \ 2 = ^data to predict 
the correct A; then the function is again minimized with the new A, and so on, until the value of the 
X 2 falls within the proper range. The minimization is performed using the Polak-Ribiere conjugate 
gradient method (Press et al. 1992). The expressions for the gradients of % 2 , H±, and H2 are given 
in Appendix B. 



s 2 + rl(sm 2 6 j + ^cos 
1 9b 



(12) 



2.3. Velocity Calculations 



To examine the dynamical effects of the projection degeneracies we assume a constant mass- 
to-light ratio, axisymmetric, two-integral dynamical model (e.g. Binney & Tremaine 1987, van der 
Marel 1991). Following Binney et al. (1990), we calculate the potential <E> from the mass density 
p = Tqu by first expanding the density in Legendre polynomials, 



Pl (r') = T / i/(r',0)Pi(cos0)siii0d0 
Jo 

and then by finding the potential produced by the Legendre expansion of the density 

dr' 



$(r,0) 



r 1 r r r°° dr 

-2vrG£PKcos0) [-^ j o Pl (ry^dr' + r l l P^r')-^ 



The solutions to the Jeans equations for the two-integral model are 



va 



(13) 



(14) 



(15) 
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and 



( „ ?)(fl ,z) = ^ + JJ _ + _J_2, 



(16) 



where the galaxy is assumed axisymmetric, steady-state ((^z) = = 0), and isotropic in the 
meridional direction (or = oz = c) (see Binney &; Tremaine 1987 §4.2). The projected line-of-sight 
velocity dispersion is 



dz 



— OO 

oo 



ua 2 (R, Z)(cos 2 i + sin 2 0sin 2 i) + v{v^) cos 2 0sin 2 i 

<9<£» d(vo 2 ~ 



/OO /* OG 

vcj 2 (R,Z)dz + s\n 2 i / 
-oo J —oo 



COS 



dz. 



(17) 
(18) 



We numerically integrate equations (15) and (16) to find the velocity components a 2 and (v^) in 
the galaxy's meridional plane, and then we numerically integrate equation (17) to find (vf } (we do 
not separate out its components {vf os ) = af os + (fi os ) 2 )- The accuracy of the numerical integration 
algorithm was verified with the analytic solution of Satoh (1980; eqs. [8] and [13]) and other more 
ad hoc analytic solutions. 



3. TEST PROBLEMS 

We test the method in four stages. First we confirm that if we bias the solution toward the true 
density, the method converges to the correct result. Next we demonstrate that the deprojection is 
degenerate by examining the range of solutions we can produce by altering the bias function. To 
show that the degeneracy is not a numerical artifact, we next investigate the effects of the numerical 
resolution, observational errors, and numerical errors. Finally we examine the dynamical effects of 
the degeneracy on the two-integral models. 



3.1. Deprojection with Correct Bias 

Before examining real galaxies, we studied the deprojection algorithm and the effects of de- 
generacy with artificial images of known density distributions. We used the axisymmetric, oblate 
(qo < 1) density distribution, 

/ R 2 z 2 l\ ~ a ° /2 
v {R,Z) = v Q [l + -p + -p- ? ) , (19) 

which has the analytically calculable surface brightness, 

^l/2-r(a {) -l\ / 2 2 i \ ( 1 -"o)/2 

J(s,y)= ,^r } , Wo [l+ X - + y -^i 1 H • (20) 

r (^) ,/fl 2 sin 2 i + cos 2 i V s o s 2 q 2 sm 2 i + cos 2 iJ 
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In the test image, the surface brightness points l[ m = I(wi,ip m ) are logarithmically-spaced in radius 
w (46 points from zu/so = 0.53 to w/so = 38.2), and equally-spaced in angle ip (7 points from 
iP = to tp = vr/2). 

The errors and uncertainties in the image of a bright galaxy are dominated by systematic effects 
rather than by photon counting noise. Imperfect flat-fielding produces a fractional error which varies 
over large scales in an image, and the sky background level is uncertain to within some constant 
value. For a set of image points J^, we adopted an error model with o\ m = gqI^ + I^. To simulate 
a real image, we added artificial errors to our "data". A constant background offset, lb, simulated 
an error in determining the sky background level, and a Gaussian distribution of fractional errors 
(with rms amplitude ao, and Fourier width corresponding to fluctuation wavelengths of about 5% 
the full-frame width of the image) was added to simulate flat-fielding errors. These model errors 
qualitatively resemble the residuals from fitting real images. 

We deprojected simulated images for various values of the model parameters (ao, So, qo, i), the 
biasing parameters (a^, Sb, %), and error levels (ao, 1^). The results of some of these deprojections 
are shown in Figure 1, where the image has the analytic form specified by equation (20), with 
cto = 3.0, so = 1-7, qo = 0.6, ao = 0.01 (corresponding to a 0.01 magnitude error), 1^ ~ 0.0004Ij^ ax 
(sky level underestimated), and i = 30° (recall that i = 0° is pole-on), with simulated errors 
added as described above. The first column of Figure 1 shows the case in which we set the bias 
function equal to the "correct" density: h = (i-e., ctb = ao, Sb = s o, % = qo)- in this 

test case, the final solution almost exactly matches the input model, with rms fractional errors of: 
((u — fcorr)/ ^con^rms = 0.009. Similar deprojections were performed over a wide range of model 
parameters, with similar results: when given the proper bias functions, the algorithm reconstructed 
the original density function. Large changes in the artificial error parameters had negligible effects. 

3.2. The Existence of Degeneracies 

Our real interest, however, is how the results can differ from the input density model given 
the degeneracies in the projection operator. Figure 1 also shows three other deprojections of the 
same image, where each column shows a solution with a different bias function h. We delib- 
erately chose extreme biases and a low inclination (i = 30°) to highlight the range of possible 
solutions. In the second column, the bias function is more elliptical than the "correct" solu- 
tion (ab = 3.0, Sb = 1.7, % = 0.3); in the third column it decreases more rapidly with radius 
(«b = 3.5, Sb = 2.0, % = 0.6); and in the fourth column we have added an exponential disk 
(h^ 1 = f corr + fed exp(— i?/i?b)sech(Z/Zb), with = 0.25, Rb/so = 5.2, and Z^/so = 0.9). Each 
model is converged until it is a statistically acceptable fit to the data (x 2 — ^duta)- The rms 
fractional differences in the density for the solutions are {(v — f C orr)/^corr)rms = 0.58, 0.35, and 
0.22, for columns 2, 3, and 4, respectively, demonstrating that the konus densities can have large 
amplitudes. These densities qualitatively resemble the disks examined by Rix &: White (1990) and 
the analytic solutions of Gerhard & Binney (1996) and Kochanek & Rybicki (1996), without the 
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discontinuity problems exhibited by those solutions at large radii. Although some of the solutions 
show systematic deviations in the outer regions of the galaxy at the two-sigma level, and the density 
structure of the disk-biased solution has some unphysical features, it should be remembered that 
we deliberately chose the bias functions to find extreme examples. 

3.3. Searching for "True" Konus Densities 

We next checked the effects on the solutions of varying the numerical resolution. A true 
projection degeneracy should exist independently of the sampling of either the image or the density 
distribution, except in the limit that the data overconstrains the density. Our standard N r x N a 
density grid has N r = 100 radial and N a = 25 angular zones. There are numerical errors in the 
projection, and if we compute the errors for the test problem in §3.1, the contributions to xV^data 
from numerical errors are ~ 0.003, 0.01, and 0.07 for i = 30°, 60°, and 90° respectively, given our 
standard error model. Although these numerical errors are much smaller than the observational 
errors (xV-^data — 1), we must examine whether the magnitude of the degeneracies is influenced 
by the resolution. 

We performed a series of tests comparing the numerical projection of the test model to the 
analytic projection of the test model, gradually reducing the density zone resolution. Reduced 
resolution increases the numerical projection errors and reduces the range of possible solutions 
because the number of degrees of freedom in the density model (N r x N a ) approaches the number 
of data points (A^ata = 322). The inclination was fixed at i = 30°. For zone resolutions of 
N r xN a = 2500, 1600, 400, 196, and 100, the resulting numerical errors were xL/^data ^ 0.003, 
0.008, 0.2, 1.2, and 6.0, respectively. The x 2 from numerical errors became unacceptable only in 
the two lowest resolution simulations, where the density distribution was overconstrained (A zoncs = 
N r x N a < A data ). 

Next we deprojected the test model (qo = 0.6) using substantially rounder (<fo = 0.9) or 
flatter (% = 0.3) bias densities. As we reduced the resolution, the allowed degeneracy towards 
low-ellipticity models was little affected, but the degeneracy towards high-ellipticity models was 
substantially reduced (see Figure 2). The amplitude of the degeneracy did not change between 
our standard resolution and the next lowest resolution. The inability to produce high ellipticity 
models at low numerical resolutions is due to the decreasing angular resolution. Note, however, 
that the degeneracies still exist even when the data overconstrains the density, a clear sign that the 
degeneracies are not due to the numerical resolution of the calculation. 

A true degeneracy should also exist independently of the observational errors. The amplitude 
of the degeneracy may increase with the amount of noise in the observations, but it should not 
vanish in the limit of no noise. We made a sequence of deprojections of a single image while 
gradually reducing the observational error levels, cjq and 1^, until the numerical errors become 
significant (Xerr/-^data ~ 0.1); lb was reduced in proportion to ctq. Figure 3 shows the deprojections 
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of the same image used in Figure 1 (ao = 3.0, so = 1.7, go = 0.6), biased to either % = 0.9 and 
% = 0.3 as a function of the inclination (90°, 60°, or 30°), with varying errors. The permitted 
range of the solutions depends strongly on the inclination and weakly on the errors. At i = 90° 
the solution is unique, with any uncertainty in the axis ratio due to the noise; the rms fractional 
density variation between the extreme solutions ((ft, = 0.9, 0.3) at "normal" noise levels (ao = 0.01) 
is ((^0.3 — t / o.9)/ l/ o.9)rms = 0.26. At i = 60° the axis ratio can be biased upwards and downwards by 
Aq ~ 0.1 about the true axis ratio; for do = 0.01, the allowed rms variation is ((^0.3— i / o.9)/ z/ o.9)rms = 
0.53. At i = 30° the axis ratio can be biased by Aq ~ 0.2 about the true axis ratio. For do = 0.01, 
the allowed rms variation is ((^0.3 — 1/ o.9)/ t 'o.9)rms = 0.82; even when the errors are at the level of 
0.002 magnitudes, the allowed rms variation is ((^0.3 — ' y o.9)/ t 'o.9)rms = 0.67. 

In all previous tests, we fit our projections to the analytic surface brightness model, so numer- 
ical projection errors limited how far we could reduce the noise. For our final test we "removed" 
these errors by using the numerical projection of the density model as the image, so that as the 
model observational noise approaches zero there is a numerical solution with xV-^data = 0. Fig- 
ure 4 shows a sequence of deprojections of the same qo = 0.6 model, with i = 30°, in which the 
solution was biased to q^ = 0.3 and % = 0.9, and the stated error level was gradually decreased 
from do = 0.01 to do = 10 -4 (no artificial noise was added); for speed, a coarser density grid was 
used (N r x N a = 52 x 13). The allowed ellipticity range in the density decreased as ao decreased, 
but at a much slower rate than the residuals in the image (Figure 5). For a true konus solution 
we should see the rms variation of the density become constant and finite at zero noise, but the 
presence of any noise always allows the variations to be larger. Gerhard & Binney (1996) refer to 
such additional degeneracies as "truncated konus densities", whose projections are never exactly 
zero but have large rms density variations for very small projected surface densities. The continued 
presence of a coherent 46% rms density variation, even in the limit that the errors are one hundred 
times smaller than typical observational data, means that for all practical purposes we have found 
true degeneracies of the projection operator rather than any numerical effect. 



The existence of a deprojection degeneracy is only of academic importance unless it has dynam- 
ical consequences. The fundamental question is whether the dynamical uncertainties introduced by 
the deprojection uncertainties can alter the conclusions of constant mass-to-light ratio, axisymmet- 
ric, two-integral dynamical models (e.g. van der Marel 1991). For a fixed mass-to-light ratio, we 
calculated the velocities (as described in §2.3) in the meridional plane, a 2 and {v^}, and the pro- 
jected mean square line-of-sight velocity (tf os ) for solutions to the analytic galaxy image (qo = 0.6); 
these velocity profiles are shown in Figure 6 as a function of inclination and bias function. The total 
luminosity (and thus the total mass) of the solutions is fixed, and the average radial distribution of 



the luminosity cannot vary a great deal, so the total meridional plane velocities y 2cr 2 + (v^) (and 

thus the kinetic energy) remain essentially constant, with variations typically ^20% at 30° and 



3.4. 



Dynamical Consequences 
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£3% at 90°; this is a consequence of the virial theorem. However, the velocity anisotropy (v^) 1 / 2 /a 
varies considerably, with larger variations at lower inclinations (£60% at 30° vs. £25% at 90°). 
Projection of these velocities results in weak variations in the line-of-sight velocities, taking the 
form of alterations in the ratio between the major and minor axis profiles. Typical variations in 
(wgj 1 / 2 at % = 90° are £10%, at % = 60° are £20%, and at i = 30° are £30%. Such variations are 
roughly comparable to typical measurement errors in line-of-sight velocities. 

4. REAL EXAMPLES: NGC 7619 AND NGC 1439 

For our final experiment we selected for deprojection two galaxies from van der Mar el (1991), 
NGC 7619 and NGC 1439; these galaxies were chosen for their small isophote twists. The photo- 
metric profiles were taken from Franx, Illingworth, & Heckman (1989a, hereafter FIHa), where the 
R-band photometry is parametrized by an "intermediate" axis (m = y/ab) profile for each galaxy, 
along with the ellipticity, isophote twist, and higher-order isophote corrections at each radial point. 
We converted the data into a series of radial profiles and their errors at 7 evenly spaced angles 
beginning on the major axis and ending on the minor axis. We incorporated the isophote twist 
and cos 4-0 terms into our data, but none of the other higher order terms (strictly speaking, the 
isophote twist, which is an indicator of triaxiality, violates our density symmetry requirements, but 
in practice its presence was of little consequence for the final solutions). We made no corrections for 
the seeing, so the density profiles inside ~ 5 arcseconds will be unreliable. The photometric uncer- 
tainties are dominated by flat-fielding errors (at the 1% level) except at the outer radii, where there 
is a sky level uncertainty of l%-3%. Both galaxies were deprojected for a variety of inclinations 
and biasing parameters. 

The kinematic data were taken from Jedrzejewski & Schechter (1989, hereafter JS) and from 
Franx, Illingworth, & Heckman (1989b, hereafter FIHb), where the velocity dispersion and rotation 
velocity are given along the major and minor axes. The line-of-sight velocities were then derived 
from the density solution as in §3.4, and fit to the kinematic data (by minimizing the \ 2 between 
them) to find the -R-band mass-to-light ratio Yr; to minimize the effects of seeing on the results, 
and to allow for direct comparison with van der Marel (1991), any velocity data point inside 4 
arcsec is not used in the fit. Several of these solutions are shown in Figures 7 and 8. 

NGC 7619 is an E2/E3 elliptical with a small (< 13°) isophote twist; van der Marel (1991) 
was not able to fit an acceptable axisymmetric two-integral dynamical model for this galaxy. The 
projected ellipticity varies from 0.16 to 0.28, with small (< 0.7%) higher-order corrections. As is 
evident from the residuals in Figure 7, there is a radial "kink" in the outer regions of the galaxy 
(r ~ 50 arcsec) which cannot be well fit by a single power-law. The deprojections can be made 
modestly more or less elliptical, with some solutions showing boxy and "SO-like" structures. The 
range of the solutions is not as large as found for the test galaxies (compare Figure 7 to Figure 3) , 
probably because much of the "smoothness" was taken up in trying to fit to a single radial power- 
law. The rms fractional density variations between the extreme solutions (% = 0.9 and % = 0.5) 
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for inclinations of 90°, 60°, and 40° are ((fo.5 — ^o.9)/^o.9)rms = 0.07, 0.15, and 0.34, respectively 
(i = 40° was used because a convergence difficulty caused minimization time to be prohibitive for 
i = 30°). For i = 90°, the solutions (with different biases %) had iZ-band mass-to- light ratios 
of (T/T ) fl ~ (3.9±0.2)/i 50 and x 2 /N ~ 2.8 for the kinematic fit, where H = 50/i 50 km s" 1 
Mpc -1 ; the error bars represent the range of solutions with A\ 2 = ±4 (where the poor fits of these 
solutions were renormalized to set x 2 = N). For i = 60°, the range of solutions had (T/T )# ~ 
(3.8-3.9±0.3)/i 50 and x 2 /N ~ 3.2-4.3. For i = 40°, the solutions had (T/T©) fl ~ (4.1±0.3)- 
(4.1±0.4)/i5o and x 2 /N ~ 3.4-5.8. Given the uncertainties and the fact that neither we nor van der 
Marel found an acceptable (x 2 /N ~ 1) solution for the velocities, our mass-to-light ratios are 
consistent with van der Marel's (T/T Q ) R ~ (4.0 ± 0.1)h 50 at 90° and (T/T ) R ~ (4.1 ± 0.1)h 50 at 
60°. 

The second galaxy, NGC 1439, is an El elliptical with a small twist angle < 9° and an indication 
of a disklike distortion in the inner parts (FIHa). The projected ellipticity varies from 0.07 to 0.11, 
with higher-order corrections < 0.8%. It has a counter-rotating core (FIHb). As can be seen in the 
residuals in Figure 8, the projected galaxy has a strong variation in ellipticity with radius which 
is not well fit by a constant ellipticity. Solutions biased toward high ellipticities again develop 
disk-like structures that sometimes show a feature along a line at the inclination angle from the 
symmetry axis. The rms fractional density variations between the extreme solutions (% = 1-0, 
q h = 0.6) for inclinations of 90°, 60°, and 30° are ((u . 6 - ^i. )M.o)rms = 0.05, 0.11, and 0.33, 
respectively. For i = 90°, the solutions had (T/T Q ) R ~ (2.3±0.3)/i 50 and x 2 /N ~ 1.7. For i = 60°, 
the solutions had (T/T ) fl ~ (2.3±0.3)/i 50 and x 2 /N ~ 1.7-1.8. For i = 30°, the solutions had 
(T/T )ij ~ (2.4-2. 5±0.3)/i5o and x 2 /N ~ 1.6-1.8. Given the uncertainties, our mass-to-light ratios 
are consistent with van der Marel's (T/T )# ~ (2.2 ± 0.1)/i5o- 

Both galaxies show similar variations of the line-of-sight velocities. The velocities in the merid- 
ional plane show large variations (^20% for NGC 1439 at i = 30°), creating differences in the 
projected velocities that are smaller (^7%) than typical measurement errors. The same geometric 
effects which create the konus densities appear to also create "konus velocities" whose large am- 
plitudes in the meridional plane practically vanish in projection; thus, considerable refinement of 
spectroscopic techniques would be needed to rule out any deprojection degeneracies with velocity 
measurements. It is possible that higher-order velocity moments are not as strongly affected by the 
konus degeneracy, but an analysis of this question was outside the scope of this paper. Note that 
none of the solutions actually fits the velocity data well. While this does not alter our conclusions 
about the degeneracy of a real deprojection solution, an accurate model would need to account for 
seeing effects, a varying mass-to-light ratio, anisotropics, and triaxiality, in order to acceptably fit 
the data. Note also that only two-integral models were used, and it is possible that three-integral 
axisymmetric models would show larger kinematic variation in projection. 
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5. CONCLUSIONS 

The deprojection of an axisymmetric galaxy is uniquely specified only if i = 90° and the 
symmetry axis is in the plane of the sky. At all other inclinations there is a gradually increasing 
degeneracy in the projection operator corresponding to an unconstrained "cone of ignorance" with 
opening angle 90° — i in the Fourier transform of the density (Rybicki 1987). Recent analytic studies 
by Gerhard & Binney (1995) and Kochanek & Rybicki (1996) have found simple density functions, 
called konus densities, that are invisible in projection because their Fourier transforms are non-zero 
only in the cone of ignorance. 

We have developed a new deprojection method based on linear regularization and explored the 
effects of the deprojection degeneracy on the inferred structure and dynamics of axisymmetric ellip- 
tical galaxies. The advantages of our approach over earlier methods are that it is non-parametric, 
that it performs a well-defined statistical fit to the surface brightness data, that it strictly enforces 
the positivity and monotonicity of the solution, and that it allows us to explore the degeneracies 
of the projection operator. The standard method of Lucy (1974) is non-parametric but does not 
have a well-defined convergence criterion, and functional fitting methods such as Palmer's (1994) 
and Bendinelli's (1991) depend on parametric forms. No previous method has been able to explore 
the degeneracies of the projection operator, or to impose monotonicity or any analytic requirement 
on the models. 

We find that axisymmetric galaxies have large deprojection uncertainties at modest inclinations 
even when we are restricted to positive definite, monotonic density distributions. The uncertain- 
ties are not due to numerical projection errors, insufficient grid resolution, or observational noise, 
although increasing the noise in the observations increases the uncertainty. Even when the obser- 
vational error in the surface brightness points approaches 10 -4 mag, it is possible to have density 
distributions fitting the data with rms fractional variations of 46% for i = 30°. The differences 
between the model densities are the konus densities, and they resemble the analytic solutions found 
by Gerhard & Binney (1995) and Kochanek & Rybicki (1996). If the bias function used to produce 
variations in the model density is not too spherical or too elliptical compared to the true density, 
the resulting model density looks reasonable. Solutions biased toward very high ellipticities show 
strong "disk-like" structures with a feature at angle i from the equatorial plane, and could be 
rejected as physical inversions. Because of the noise and konus degeneracy, the solutions are quite 
sensitive to the choice of the bias function; we chose simple functions which do not reproduce well 
some of the more complicated features in the data (e.g. NGC 7619's radial kink), but one could 
easily implement a more "accurate" bias by using Lucy's method or a parameterized approximate 
deprojection method to arrive at an initial bias function. 

We have evaluated the dynamical variations allowed by the deprojection uncertainties in the 
constant mass-to-light ratio, axisymmetric, two-integral dynamical model (Binney et al. 1990; van 
der Marel 1991). Although the velocities in the meridional plane can have large variations, the 
variations in the projected mean square velocities are modest for all inclinations when compared 
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to typical measurement errors. Given our structural and dynamical assumptions, current velocity 
measurements are not helpful in reducing the deprojection degeneracy — the konus densities are 
associated with what we might term "konus velocities" . We infer mass-to-light ratios for NGC 1439 
and NGC 7619 which are comparable to those of van der Marel (1991), but with larger uncertainties. 
None of the constant mass-to-light models fits the data well, with typical xV^dof — 1-5-3. 

The implications of our results for more complicated models are unclear. Statler (1994a, 
1994b) and Statler & Pry (1994) have used dynamical models to tightly constrain the deprojection 
of a triaxial galaxy, but they have assumed axis-ratios constant with radius, severely limiting the 
generality of the solution. Merritt (1996) has developed a technique for deriving a unique two- 
integral distribution function from surface brightness and velocity moment measurements, but only 
at the non-degenerate edge-on inclination; presumably, the large uncertainty in the density due to 
the degeneracy must also affect the inferred distribution function. 



ACKNOWLEDGMENTS 
We thank George Rybicki, Dan Fabricant, and Roeland van der Marel for their helpful comments. 



-15- 



REFERENCES 

Arnold, R., Robijn, F. H. A., k de Zeeuw, P. T. 1995, in Fresh Views of Elliptical Galaxies, ed. A. 
Buzzoni, A. Renzini, k A. Serrano (San Francisco: ASP), 23 

Bendinelli, O. 1991, ApJ, 366, 599 

Bertola, F., Bettoni, D., Danziger, J., Sadler, E., Sparke, L., k de Zeeuw, T. 1991, ApJ, 373, 369 
Binney, J. J., Davies, R. L., k Illingworth, G. D. 1990, ApJ, 361, 78 
Binney, J., k Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton Univ. Press) 
Carollo, C. M., k Danziger, I. J. 1994, MNRAS, 270, 743 

Carollo, C. M., de Zeeuw, P. T., van der Marel, R. P., Danziger, I. J., k Qian, E. E. 1995, ApJ, 
441, L25 

Dehnen, W. 1995, MNRAS, 274, 919 
Dehnen, W., k Gerhard, O. E. 1993, MNRAS, 261, 311 
Dehnen, W., k Gerhard, O. E. 1994, MNRAS, 268, 1019 
Fabbiano, G. 1989, ARA&A, 27, 87 

Fabbiano, G. 1995, in Fresh Views of Elliptical Galaxies, ed. A. Buzzoni, A. Renzini, k A. Serrano 
(San Francisco: ASP), 103 

Franx, M., Illingworth, G., k de Zeeuw, T. 1991, ApJ, 383, 112 

Franx, M., Illingworth, G., k Heckman, T. 1989a, AJ, 98, 538 (FIHa) 

Franx, M., Illingworth, G., k Heckman, T. 1989b, ApJ, 344, 613 (FIHb) 

Franx, M., van Gorkom, J. H., k de Zeeuw, T. 1994, ApJ, 436, 642 

Gerhard, O. 1991, MNRAS, 250, 812 

Gerhard, O., k Binney, J. 1996, MNRAS, in press 

Hunter, C., k Qian. E. 1993, MNRAS, 262, 401 

Jedrzejewski, R. k Schechter, P. L. 1989, AJ, 98, 147 

Kochanek, C. S. 1995, ApJ, 445, 554 

Kochanek, C. S. 1996, ApJ, in press 

Kochanek, C. S., k Rybicki, G. B. 1996, MNRAS, in press 

Lucy, B. 1974, AJ, 79, 745 

Merritt, D. 1996, AJ, in press 

Maoz, D., k Rix, H.-W. 1993, ApJ, 416, 425 

Newton, A. J., k Binney, J. 1984, MNRAS, 210, 711 

Palmer, P. L. 1994, MNRAS, 266, 697 



-16- 



Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes in 
C, 2nd Ed. (Cambridge: Cambridge Univ. Press) 

Qian, E. E., de Zeeuw, P. T., van der Marel, R. P., & Hunter, C. 1995, MNRAS, 274, 602 

Rix, H.-W., & White, S. D. M. 1990, ApJ, 362, 52 

Rybicki, G. B. 1987, in Structure and Dynamics of Elliptical Galaxies, ed. P. T. de Zeeuw (Dor- 
drecht: Kluwer), 397 

Saglia, R. P., et al. 1993, ApJ, 403, 567 

Satoh, C. 1980, PASJ, 32, 41 

Schweizer, F., Whitmore, B. C, & Rubin, V. C. 1983, AJ, 88, 909 
Statler, T. S. 1994a, ApJ, 425, 458 
Statler, T. S. 1994b, ApJ, 425, 500 

Statler, T. S. 1995, in Fresh Views of Elliptical Galaxies, ed. A. Buzzoni, A. Renzini, & A. Serrano 
(San Francisco: ASP), 27 

Statler, T. S., & Fry, A. N. 1994, ApJ, 425, 481 

van der Marel, R. P., Binney, J., & Davies, R. L. 1990, MNRAS, 245, 582 
van der Marel, R. P. 1991, MNRAS, 253, 710 

van der Marel, R. P., Evans, N. W., Rix, H.-W., White, S. D. M., & de Zeeuw, T. 1994, MNRAS, 
271, 99 

Whitmore, B. C, McElroy, D. B., & Schweizer, F. 1987, ApJ, 314, 439 



This preprint was prepared with the AAS IAT^X macros v4.0. 



-17- 



APPENDICES 



A. First-Order Projection Scheme 

From the interpolated density of equation (6), we find its contribution to the projection: 



dljkim = / Vjk + — —r k jdz- J —r(z)dz 

+ ^±i_S rt ) ( , _ fl) - h!±i^2i I" rtz)iz , (Ai) 

r k - r k+1 J r k - r fc+1 J Zl 



Since r 2 = w 2 + z 2 , 



Jrdz = J Z \w 2 + z 2 ) 1 / 2 dz 

2 

= Z -^Jw 2 + z 2 + ^- In (z + Vw 2 + z 2 ) 



zr w 2 , , ... 
T + — ln(z + r), (A2) 



which leads to the first-order approximation in equation (7). 



B. Gradient Expressions 

The derivatives of the \ 2 statistic are: 
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The derivatives of the radial and angular smoothness parameters H\ and H 2 are: 
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CAPTIONS 



Fig. 1. — Effects of the bias function on deprojections of an artificial galaxy image. Each column 
represents a different choice for the bias function. The top row shows the normalized residuals of 
the solution, (7 — I*) /a, where the horizontal dotted lines are the one-standard-deviation errors. 
The dashed lines are the residuals of the projection of the "correct" density (i.e. they show the 
added "noise"), and the curves are offset for visibility with the major axis on top. The second row 
shows the deprojected density profiles on fixed azimuths. The third row shows contour plots of the 
density solution (solid contours) and the "correct" density (dotted contours). There are 26 contour 
levels shown, logarithmically spaced from v = 0.25 to v = 2.5 x 10~ 6 . The bottom row shows 
contour plots of the konus density, or the difference between the final solution and the "correct" 
density; dotted contours represent negative values. There are 42 contour levels shown (positive and 
negative), logarithmically spaced from v = 1CT 2 to v = 1CT 6 . 

Fig. 2. — The effects of numerical resolution on deprojections of artificial galaxy images. The 
inclination is fixed at i = 30°. The major-to-minor axis ratio is plotted for each solution as a 
function of radius. The "correct" solution is the central dashed line with qo = 0.6; the upper 
and lower sets of curves are the solutions when biased toward % = 0.9 and q^ = 0.3. The solid, 
dashed, dotted, and dot-dashed curves show the solutions when the density zone resolution is 
N r x N a = 100 x 25, 80 x 20, 40 x 10, 28 x 7, and 20 x 5, respectively. Note that the image 
resolution is N x x N y = 46 x 7, so the density is overconstrained by the data in the last two cases. 

Fig. 3. — Deprojections of an artificial galaxy image as a function of noise level and inclination. 
Each column represents a different assumed inclination i, with i = 90°, 60°, and 30°, from left to 
right. In the two top rows, the solid, dotted, and dot-dashed lines show solutions with decreasing 
fractional noise aQ. The highest noise level is set to a realistic value; the lowest is the level at which 
numerical errors become significant. The top row shows the density axis ratio, where the central 
dashed line is the deprojection using the "correct" bias. The middle row shows the density solution 
angular profile at a fixed radius (r/sQ = 1.6). The bottom row shows contours of the solution 
density. The dotted and dashed lines show solutions at the highest (normal) noise level (% = 0.9 
and % = 0.3, respectively); the solid contours show the "correct" solution, qo = 0.6. There are 26 
contour levels shown, logarithmically spaced from v = 0.25 to v = 2.5 x 10~ 6 . 

Fig. 4. — Deprojection of an artificial galaxy image with numerical errors removed (see text), as 
a function of the noise level uo- The inclination is fixed at i = 30°. The axis ratio is plotted for 
each solution. The "correct" solution is the central dotted line with c/o = 0.6; the upper and lower 
sets of curves are the solutions when biased toward % = 0.9 and % = 0.3. The solid curves show 
a sequence of solutions with noise level uo lowered. 

Fig. 5. — Effects of observational errors on deprojection. The rms fractional variation in the density 
solution and the rms fractional error in its projected image are shown as a function of error level. 
The dotted line solutions were biased toward c/b = 0.3, and the solid line solutions were biased 



-20- 



toward = 0.9. 

Fig. 6. — Velocity dispersions for an analytic test galaxy (qo = 0.6), assuming a constant mass- 
to-light ratio. The left column shows the velocities in the meridional plane, along the major axis 
(the variations along the minor axis were inconsequential). The right column shows the projected 
velocities along the major (upper profiles) and minor (lower profiles) axes. The top box shows 
solutions for i = 90°, the middle for 60°, and the bottom for 30°. At each inclination, solutions are 
shown with biases of % = 0.9, 0.6, and 0.3. The vertical axis has arbitrary units. 

Fig. 7. — Deprojection solutions and resultant projected velocities for NGC 7619, for varying 
inclination and bias function. Each column represents a different assumed inclination i. The top 
row shows the normalized residuals of the solution: (I — I*) /a, where the horizontal dotted lines 
are the one-standard-deviation errors. The middle row shows the rms line-of-sight velocity profile 
for the major (upper) and minor (lower) axes. Velocity data (JS; FIHb) are superimposed, where 
the crosses are major axis points, the solid squares are minor axis points, and the error-bars are 
one-<7 errors. The bottom row shows contours of the solution density. There are 28 contour levels, 
from v = 0.25 to v = 10 -6 . The solid lines indicate a solution with a % = 0.7 bias function; the 
dotted lines, % = 0.9; the dashed lines, q^ = 0.5. The vertical dot-dash line indicates the cut-off 
radius (4 arc sec) below which the velocity data are not used. 

Fig. 8. — Deprojection solutions and resultant projected velocities for NGC 1439, for varying 
inclination and bias function. Velocity data are from FIHb. Labels are the same as in Fig. 7. 
There are 24 contour levels, from v = 0.25 to v = 6.3 x 10~ 6 . The solid lines indicate a solution 
with a <7b = 0.8 bias function; the dotted lines, % = 0.99; the dashed lines, % = 0.6. 
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